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We study numerically the distribution of scattering phases V($) and of Wigner delay times V(jw) 
for the power-law banded random matrix (PBRM) model at criticality with one channel attached to 
it. We find that - P($) is insensitive to the position of the channel and undergoes a transition towards 
uniformity as the bandwidth b of the PBRM model increases. The inverse moments of Wigner delay 
times scale as (r^> ~ L~ qD "+ 1 , where D q are the multifractal dimensions of the eigenfunctions of 
the corresponding closed system and L is the system size. The latter scaling law is sensitive to the 
position of the channel. 
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I. INTRODUCTION 

The transport properties of random media have been a 
subject of intensive research activity for more than fifty 
yeara 1 ! 2 !?! 4 !^!?!?!?! 1 ?! 11 ! 12 ! 1 ?. The most fascinating as- 
pect of this study is the appearance (in high-dimensions) 
of a metal-insulator phase transition (MIT) as an ex- 
ternal parameter changes. In the metallic phase, the 
eigenstates are extende d 3 i 4 i 7 i 8 and the statistical prop- 
erties of the spectrum are described well by random ma- 
trix theory 2,3,4 *^. In the localized phase, the levels be- 
come uncorrelated leading to a Poissonian level spacing 
distribution^, and the eigenfunctions are exponentially 
localized 1 ! 2 ! 3 ! 4 ! 7 . 

The MIT where the phase transition from localized 
to extended states occurs, is characterized by remark- 
ably rich critical propertiesiSiSiSiiSiiiii^. In particular, 
the eigenfunctions show strong fluctuations on all length 
scales and represent multifractal distributions 7 ' 8,9 . This 
is one of the most important features at MIT and has con- 
siderably helped our understanding of different phenom- 
ena in mesoscopic systems 2 . The multifractal structure of 
the eigenfunctions is usually quantified by studying the 
size dependence of the so-called participation numbers 
(PN) 



Afq = 



\ip{r)\ 2q dr I cx 



(1) 



where L is the linear size of the system and D q are 
the multifractal dimensions of the eigenfunction ^(r). 
Among all the dimensions, the correlation dimension D2 
plays the most prominent role. The corresponding PN 
is roughly equal to the number of non-zero eigenfunction 
components, and therefore is a good and widely accepted 
measure of the extension of the states. At the same time, 
D2 manifest itself in a variety of other physical observ- 
ables. As examples we mention the conductance dis- 
tribution in metals^i^, the statistical properties of the 
spectrum 1 ^, the anomalous spreading of a wave-packet, 
and the spatial dispersion of the diffusion coefficienlii 3 .. 

It is therefore highly desirable to have the possibility 
of measuring the multifractal dimensions D q , and specifi- 
cally D2, directly from an experiment. Unfortunately, to 



our knowledge, there is no experimental study of eigen- 
function fluctuations at the transition regime—, although 
in recent years a MIT was experimentally accessible in 
microwave 1 ^ and optical wavei 7 - setups. The main diffi- 
culty is the lack of eigenfunction accessibility in all these 
experiments. On the other hand, new techniques were 
developed that allow accessibility to various quantities 
associated with the scattering process. Among them is 
the Wigner delay timei 8 ^* 2 ^ which captures the time- 
dependent aspects of quantum scattering. It can be in- 
terpreted as the typical time an almost monochromatic 
wave packet remains in the interaction region. It is re- 
lated to the energy derivative of the total phase shift 
$(_E) = — imdet S(E) of the scattering matrix S(E), i.e. 
tw(E) = jj rf ^f ^ , where M is the number of channels. 

These experimental efforts have promoted the study of 
Wigner delay times to quantities of interest in their own 
right and attracted a lot of theoretical work. For chaotic 
systems many results are known concerning the distribu- 
tion of Wigner delay times V(tw)^- Recently, the inter- 
est has extended to systems showing diffusion's 3 , and 
localization 24 . At the same time, an intensive activity 
to understand V(tw) for systems at critical conditions 
was undertaken inSi 2 ^ 2 ^ 7 .. A possibility of anomalous 
scaling of inverse moments of Tyy with the system size L 
was suggested inS 3 .. Specific predictions linking the scal- 
ing exponents and the multifractal properties of eigen- 
functions of the corresponding closed system were made 
in 25,26 . While in both cases an anomalous scaling of the 
inverse moments of delay times was reported, there was 
a discrepancy in the exact power dictating these scaling 
laws. 

Precise reasons for these discrepancies are still not 
clear (see though the discussion inSi) and deserve in- 
vestigation. Here we undertake this task by studying 
numerically the distribution of Wigner delay times at 
MIT for the prototype power-law banded random ma- 
trix (PBRM) model with one channel attached to it. We 
provide clear evidence showing that the inverse moments 
of Wigner delay times (t^ 9 ) scale as 



Mi- 



ce L 



-fit) 



/(«) = qD 



9+1 



(2) 



where (.) stands for an ensemble average in agreement 



2 



with the theoretical prediction of 26 . However we find that 
this relation is extremely fragile. Namely we show that 
it holds for channels attached to a typical position in the 
sample. At the same time we perform a detail analysis 
of the distribution of phases V{$) of the scattering ma- 
trix associated with our scattering setup. We show that 
■p($) is insensitive to the position of the channel while 
it depends drastically on the band-width b which char- 
acterize the criticality. It undergoes a transition towards 
uniformity as b increases. We note that the theoretical 
analysis o& applies only for b 3> 1 while our numerical 
investigations extends over the whole &-range. 

The structure of this paper is as follows. In the next 
section we describe the scattering setup and summarize 
the properties of the corresponding closed system. Sec- 
tion III discuss the distribution of phases for various 
bandwidths, coupling strengths, and position of the chan- 
nel. In Section IV we analyze the statistical properties of 
delay times. Finally, our conclusions are given in Section 
V. 



II. SCATTERING SETUP 

The isolated sample is represented by L x L real sym- 
metric matrices whose entries are randomly drawn from 
a normal distribution with zero mean and a variance de- 
pending on the distance of the matrix element from the 
diagonal 



<(^) 2 ) 



1 



l + {\i-]\/b) 2a 



(3) 



where b and a are parameters. In a straightforward inter- 
pretation, the PBRM model describes a ID sample with 
random long-range hopping. Also, the PBRM ensemble 
arises as an effective description to a variety of physi- 
cal systems, such as the quantum Fermi accelerator, the 
scattering by a Coulomb center in an integrable billiard, 
and in billiard systems with non-analytic boundaries. 

Field-theoretical considerationaiSS^S and detail nu- 
merical investigationsiSiii verified that the model un- 
dergo a transition at a = 1 from localized states for 
a > 1 to delocalized states for a < 1. This transition 
shows all the key features of the Anderson MIT, includ- 
ing multifractality of eigenfunctions and non-trivial spec- 
tral statistics at the critical point. At the center of the 
spectral-band a theoretical estimation for the multifrac- 



tal dimensions D q gives^ 



D„ 



4&r( 9 -i/2)[V5F(«-i)r( ? -i)] 

1 - q(2Trb)- 1 



- 1 ,6<1 
,6» 1 



(4) 

where T is the Gamma function. Thus model |PJ pos- 
sesses a line of critical points b £ (0, oo), where the mul- 
tifractal dimensions D q change with b. 

We turn the isolated system to a scattering one by at- 
taching one semi-infinite single channel lead to it. The 



lead is described by a ID semi-infinite tight-binding 
hamiltonian 



#icad = >< n+ 1| + |n + 1 >< n\) . (5) 

n=l 

Using standard methods^ one can write the scattering 
matrix in the forrn^i 



S{E) = e^ {E) = l-2urW' 



1 



E — Hcs 



W, (6) 



where 7Y c ff is an effective non-hermitian Hamiltonian 
given by 



H ct i = H -inWW T . 



(7) 



Here W nm — Wo5 nno is a L x 1 vector where wq is the 
coupling strength and no is the site at which we attach 
the lead. All our calculations take place in an energy 
window close to the band center (E = 0). The delay 
time is given by 21 " 2| 



tw(E = 0) 



dE 



-21mTr(E - H eS )~ 



E=Q 



\E=0 



(8) 

We will see latter on that the scaling properties of in- 
verse moments of delay times are independent of wq while 
they depend drastically on the position hq. Specifically, 
if no is a representative site (like any site in the bulk of 
the sample) where ergodicity of wavefunctions can be as- 
sumed, then Eq J2J applies. Instead if n belongs to the 
boundary (i.e. uq — 1 or L), then Eq Q does not valid 
any more. 

In our numerical experiments bellow we used matrices 
of size varying from L = 50 up to L — 2000. For bet- 
ter statistics a considerable number of different disorder 
realizations has been used. In all cases we had at least 
100000 data for statistical processing. 



III. DISTRIBUTION OF PHASES 

We start our analysis with the distribution P($) of 
phases of the S— matrix. Analytical considerations ap- 
plied to the large b ^> 1 limiljSi led to the following ex- 
pression 



1 



27r(7 + ^Ji 1 - lcos($)) 



(9) 



where 7 = (1 + |(S)| 2 )/(1 - K^)! 2 )- Notice that in the 
limit of I (S) I = we recover the uniform distribution, 
where the corresponding setup is associated with perfect 
coupling of the sample to the leads. Eq. © was also 
derived in the framework of Random Matrix Theory^ 
using the maximum entropy ansatz. It is interesting that 
the same expression applies (see bellow) in our case as 
well. 



3 




FIG. 1: Left panels: \{S)\ as a function of wo for a channel 
attached to the first site (boundary) of the sample. Right pan- 
els: 7 3 ($) corresponding to wo where \{S)\ takes its minimum 
value (see left panels) and for various system sizes, ranging 
from L = 50 up to L = 800. In the insets of the left panels 
we show a representative V(&) for \(S)\ ~ 0.5. The dashed 
line is Eq. @. From top to bottom: b = 0.1, 1, and 10. 



In Fig. ^ we report our data for various values of b and 
various coupling strengths wq. The scattering setup cor- 
responds to one channel attached to the boundary of the 
sample. For large b values, T(<&) agrees quite good with 
the theoretical prediction of Eq. @ while for b < 1 we ob- 
serve large deviations for any wq, even for | {^(wo)) | ~ 0. 
For KS^wo))! ^ and small values of 6, two peaks ap- 
pear around <f> = it while a gap is created between them. 
As b increases further, the two peaks move closer to one 
another forming one central peak at $ = ir. In the limit 
of KS^wo))! — * 1 the distribution becomes a 5— function 
signifying a total reflection at the point where the lead is 
attached. Similar results (not reported here) are found 
for the case where the lead is attached to the bulk of the 
sample. 



IV. DISTRIBUTION OF DELAY TIMES 

In this section we analyze the distribution of Wigner 
delay times (JHJ). Recent investigations 25,26 suggested that 
inverse moments of tw scale in an anomalous way with 



respect to the sample size L. It was further suggested 
that the scaling exponent is associated with the mul- 
tifractal dimensions of the eigenfunctions of the corre- 
sponding closed system. This approach allow for an ex- 
perimental measurement of the fractality of the wave- 
functions with scattering techniques. However, these 
contributions concluded on a different scaling exponent. 
It is therefore desirable to investigate the origin of the 
discrepancy and specify the conditions under which the 
theoretical expression J2J is applicable. 

In Fig. [3 we report the delay time distribution 
^(lnr^ 1 ) for the case where the lead is attached to the 
bulk of the sample. Similar distributions (not shown) 
are found for the case where the lead is attached to the 
boundary of the sample. We report results for various 
system sizes for band widths b = 0.1, 1, and 10. The 
shape of the distributions for a specific b— value changes 
with L, leading to the conclusion that various moments of 
the distribution of inverse delay times scale in a different 
way with the system size. In the insets we concentrate 
on the scaling of the first inverse moment (t^ 1 ) ~ L~ a . 
A nice scaling is observed. Using a least square fitting 
we can extract the scaling exponent a. 

Our results are summarized in Fig. [3] where we com- 
pare the scaling exponents obtained for the two scatter- 
ing setups (channel attached at the boundary and to the 
bulk) together with the theoretical expectation which 
predicts that in the limit of b 1 the scaling exponent 
a = £>2 ■ For completeness we include also the value of D2 
numerically extracted from the scaling analysis of eigen- 
vectors of the closed systemii. We observe that Eq. (J2J 
for the (t^ 1 ) valids quite well as long as the channel 
is connected to the bulk of the sample where the ergodic 
hypothesis for the statistical properties of the cigcnstatcs 
of the corresponding closed system holds. Although the 
theoretical expression J2J was derived under the assump- 
tion of b 3> 1 we see that our numerical results for & < 1 
can also be described by Eq. relatively good. More- 
over, we find that the scaling relation J3J) persist for any 
value of (S) ^ in agreement witbJ&. On the contrary, 
for the case where a channel is attached to the boundary, 
we observe strong deviations from Eq. (J5J|. This is be- 
cause the boundary is not a representative position and 
the ergodicity hypothesis as far as the eigenfunctions is 
concerned is not any more applicable. 



Let us finally comment on the anomalous scaling be- 
havior of the other moments (t^ 9 ). A simple homoge- 
neous fractal would be completely characterized by two 
of the moments and the respective f(q) curve would be a 
straight line while a non-linearity in f(q) would signify a 
multifractal behavior. In Fig. 01a) we report our results 
for the case with a channel attached to the boundary. 
Again we see that the numerical data deviates from the 
theoretical predictions (0J for any value of b. Instead, the 
agreement is quite good for the case where the channel 
is attached to the bulk of the sample (see Fig.^fc). 
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FIG. 2: lnT^lnr^ ) as a function of lnr^ scaled to Inlr^ ) 
for a channel attached to a representative site (bulk) of the 
sample, (a) 6 = 0.1, (b) 6 = 1, and (c) b = 10. In the insets 
the scaling (t^ 1 ) oc L~ a is shown together with a linear fitting 
used to extract a. 



A a (boundary l<S>~0l) from <T W ~ > 

■ a (bulk l<S>l~0.5) from <t w '> -~y'' 
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FIG. 3: Scaling exponent a and correlation dimension of 
eigenfunctions D2 as a function of b. a is extracted from the 
scaling {t^, 1 } oc L~ a for a channel attached to the bulk and 
to the boundary of the sample. D2 calculated theoretically 29 
and numerically— from A/2 is also shown, see eq. Q. 
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V. CONCLUSIONS 

In summary we have studied the distribution of scat- 
tering phases V($>) and delay times for the power law 
banded matrix model at criticality. For small bandwidths 
b <C 1 the distribution of phases is non-uniform irrespec- 
tive of the value of | (S) | , while for large b the theoretical 
prediction o^ applies. 

We found that the multifractal nature of the eigenfunc- 
tions is reflected in the scaling of the inverse moments of 
Wigner delay times (t^ 9 ) with the system size L provided 
that the channel is attached to a representative position 
in the sample. We expect that these results will provide 
a new method for evaluating D2 (and in general D q ) in 
microwave and light wave experiments where tw can be 



FIG. 4: In f(q — 1) as a function of q for a channel attached 
(a) to the boundary and (b) to the bulk of the sample when 
\{S) \ ~ (full symbols). In (b) we also show ln/(g — 1) for 
\{S)\ ~ 0.5 (open symbols). The curves are the theoretical 
prediction of eq. @. f(q — 1) = (g — l)D q , see eq. 



extracted even in the presence of weak absorption. At 
the same time the appearance of the anomalous scaling 
of 

( T w) can be used as a criterion for detecting MIT. 
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